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Abstract 

We propose a method to describe the short-distance behavior of 
an interface fluctuating in the presence of the wedge-shaped substrate 
near the critical filling transition. Two different length scales deter- 
mined by the average height of the interface at the wedge center can 
be identified. On one length scale the one-dimensional approximation 



of Parry et al. [10| which allows to find the interfacial critical expo- 
nents is extracted from the full description. On the other scale the 
short-distance fluctuations are analyzed by the mean-field theory. 

PACS numbers : 68.45.Gd, 68.35.Rh 
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I. Introduction 



The analysis of physical systems usually involves some reductions in the 
description of the state of the systems. This is also the case for inhomo- 
geneous systems consisting of two coexisting bulk phases separated by an 
interface fluctuating in the presence of a substrate. Certain properties of 
such systems, e.g. those related to the adsorption phenomena, can be conve- 
niently described with the help of a single mesoscopic variable which is the 
distance of the interface from the substrate. 

In this paper we consider such a system. The substrate has the form of an 
infinite wedge extending along the y-direction with the opening angle 2ip, 
see Fig.l. The quasi-bulk phase adsorbed on the substrate is denoted as the 
/?-phase while the phase far above the substrate is denoted as the a-phase. 
The shape of the substrate is given hy z = \x\ cot ip and y) describes the 
distance of the interface from this substrate. Recently, it was pointed 
out [1-7] that the above system may undergo the critical transition in which 
the position of the central part of the interface (above the edge of the wedge) 
moves to infinity while the asymptotic parts of the interface corresponding 
to \x\ — > oo remain close to the substrate. This interfacial transition is called 
the filling transition to distinguish it from the wetting transition taking place 
on planar substrates ^. Thermodynamically the filling transition point 
is located at the bulk a-fi coexistence and the filling temperature (which 
depends on the wedge opening angle ip) is denoted as T^; T^, < T^, where 
is the wetting temperature on the planar substrate. 

In their recent paper Parry et al. used the transfer-matrix method to 



evaluate - among others - the values of the critical indices associated with the 
interfacial behavior near the filling transition. For this purpose the next step 
in the reduction of the description was made. The two-dimensional interface 
i{x, y) was replaced by the one-dimensional mid-point line £(?/) = £(0, y) (see 
Fig.l) for which the appropriate Hamiltonian was proposed. 
If, however, one is interested in the full two-dimensional structure of the 
fluctuating interface near the critical filling transition, i.e. also in the short- 
distance behavior which is not included in the reduced description then - 
at least in principle - one has to go beyond the mean-field analysis. In this 
paper we propose how to describe the two-dimensional interface close to the 
filling transition in the system with short-ranged forces 0. We expect that 
the geometry-dependent effects are important for short distances. Since the 
mid-point height does not vary too much on the short length-scale the idea is 
to fix this height at some arbitrarily chosen point and to assume the mean- 
field profile of the interface in the vicinity of this chosen point (along the 
x-direction) . Then one uses the mean-field approximation to describe the 
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Figure 1: The wedge geometry and the fluctuating a — (3 interface 

"relative" fluctuations around the fixed point. The two-point height distri- 
bution function for neighboring points consists of two parts: the one-point 
distribution corresponding to one of the points (or the average height of 
them) and the conditional probability distribution in the form of a Gaussian 
with position-dependent dispersion. Such quantity does not diverge at the 
filling transition and may turn out useful when some geometry-dependent ob- 
servables are considered. The mean-field description becomes then legitimate 
because by fixing the position of the interface and looking at the conditional 
distribution one forces the local fluctuations to be small and so one insists 
that the system is locally outside the critical region. 



The intcrfacial Hamiltonian in the case of a very open wedge (cot </? <^ 1) 
has the standard form [6,10] 



where i{x, y) (Fig.l) denotes the width of the adsorbed /9-like layer measured 
in the vertical direction and S is the a-P interfacial tension. a;(£) denotes the 
interfacial pinning potential corresponding to the critical wetting in the pla- 



II. The mean-field description 




(2.1) 
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nar case. For short-range forces considered in this paper it has the following 
form [6-11] 

uj{e) = -Wt exp(-£/0 + U exp(-2£/0 , (2.2) 

where ^ is the bulk correlation length (in the /3-phase), U and W are positive 
constants. (We use the convention in which the factor l/Zc^T is included 
into the Hamiltonian.) The parameter t denotes the dimensionless deviation 
from the wetting temperature for the planar substrate, i.e. t > for T < 
and t = for T = T^. £oo is the equilibrium width of the adsorbed layer on 
the planar substrate which minimizes the potential uj{i): exp(— £oo/0 = jfj'^- 
Because the wedge is very open we put sin 99 = 1 and cot ip = cos ip = a. 

The mean-field profile i{x) varies only in the x direction. It satisfies the 
Euler-Lagrange equation 

Ei"{x) = uj'{I) (2.3) 

and the boundary conditions: £(±00) = ioo, ^'(0±) = TC(- The solution of 
Eq.(2.3) is 

x{i) = ± / , (2.4) 

Ji v'2R£P7J(U)7s 

where the width of the mean-field profile at the center of the wedge io = i{0) 
satisfies uj{io) — uj{ioo) = Sa^/2. With the help of the Young equation one 
can relate uj{ioo) to the contact angle 9 on the planar substrate: —u{ioo) = 
S0^/2. From this we see that u}{io) = S(q;^ — 6^)/2 and the filling transition 
(£0 ^ 00, ioo - finite) takes place when 9(T = T^) = a. 
For small deviations 6i{x, y) = i{x, y) — £{x) from the mean field profile i{x) 
the fiuctuation Hamiltonian Hfi[Si\ = H[i + 6i] — H[i] is bilinear in 6£ 



(2.5) 



The important feature of the critical filling transition is the existence of 
the translational mode, i.e. the fiuctuation of the interface which requires 
very small energy (decreasing to at the filling point). This fiuctuation 
takes the form 5l{x) = e\i'{x)\ and the corresponding energetical cost is 
Hfi[6i] = e^aci;'(£o); it decreases to for Eq 00. 

The corresponding differential equation for correlation function G{r, r') = 
{6i{r)6i{r')) has in the mean-field approximation the following form 

[-SAr + uj"{i)] G{r, r') = 6{r - r') . (2.6) 
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However, the mean-field description fails in case of the critical filling tran- 
sition for short-ranged forces because Eq.(2.6) implies strong anisotropy of 
fluctuations of the interface. The fluctuations along the wedge diverge much 
faster than across the wedge. The latter are bounded by the geometry of 
the substrate. As shown in |]10|, the mean-field predictions are valid only for 
power-law forces of the type u!{£) ~ for p < 4. 



III. The reduction of the order parameter 



The effective way to analyze the critical fluctuations of i{x, y) near the 
filling transition point is to reduce the interfacial description by looking only 
at the mid-point height £{y) = i{0,y) [10]. In order to derive the corre- 
sponding Hamiltonian we proceed as follows: we minimize the Hamiltonian 
in Eq.(2.1) similarly as in the mean- field method but now with the constraint 
i{0,y) = i{y) imposed independently at each y |]T^|. From the corresponding 
Euler-Lagrange equation one obtains 



x{i,y) = ± 



d£i 



i V2(^(£i)-a;(£oo))/S 



(3.1) 



As the result the one-dimensional Hamiltonian Hi[i{y)] = H[i{x,y)] valid 
for configurations given in Eq.( |3.1| ) takes the form 



mm 



dy 



(3.2) 



d£i 



For short-range forces, see Eq. 
evaluated 



a;(£(7/))-cu(4 
v/2S(cu(£i)-cu(£oo))-«S 

, the above Hamiltonian can be explicitly 

s(ny))^^~-e(i-exp(-£7o) 



Ay 



e [1 - exp (-£/0]' 
+2S[(e - a)£ - 0^(1 - exp (-f/O)]" 



(3.3) 



where ^{y) = i{y) — £oo- For temperatures close to the filling transition 
one has I ^ ^ and £ ^ £oo and Eq.(3.3) reduces to the one-dimensional 
Hamiltonian proposed phenomenologically in [|l( 



H^[l{y)] ^ / dy 



e 



{e{y)Y + 2S(0 - ay{y) 



(3.4) 
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The above Hamiltonian has relatively simple structure and is easy to renor- 
malize. After introducing the rescaled variables L and Y 



Qy = (2S)-i/2(e/a - 1)^3/4^ 
it takes the form 

H,[L{Y)] = [dV 



e={2^)-'/\Q/a-l)-^/^L (3.5) 



L{Y) 



{ViY)f + LiY) 



(3.6) 



which is free from parameters. Accordingly, the critical behavior of the mean 
mid-point height {i{y)) and the correlation length C,y follow directly from the 
above rescaling: {i{y)) ~ (9 — a)~^^^ and ~ (© — a)~^^^. The values of 
the critical indices agree with those obtained in [10]. 

The one-dimensional model described by the Hamiltonian in Eq. (|3.6|) can 
be solved via the transfer- matrix method [jlO|, [Mf. However, in this method 
the presence of the factor L{Y) in front of {L'{Y)y is the source of ambiguity 
while discretizing the problem and defining the measure which is then used 
to evaluate the relevant propagator In order to avoid such problems it 
is convenient to introduce the new variable r] = 2L^^'^ /3 which "absorbs" the 
dangerous factor L{Y) in front of {L'{Y))'^. Then the Hamiltonian takes the 
form 



HMY)] = J dY [(r^'(F))V2 + (3^2)^/^] • 



(3.7) 



The corresponding propagator 



can be evaluated by solving - within the transfer matrix approach - the 
following equation 



dV_ 
dY 



(2r^2/3)^/V. 



(3.9) 



This equation must be supplemented by the appropriate boundary condition 
for rj2 = 0. The general form of such a condition 



drj2 In y 1^2=0 



at 



is similar to that found in \T^ for 2D wetting. In the present case the param- 
eter 7] is t-dependent so one expects the t-independent boundary condition 



8.3/ 2 In V 



(2S)3/4(e/a-l)3/«at. 
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For a < the edge effects become dominant and no filling is observed. Thus 
we assume a > (a~^/^ is the range of the influence of the edge effects), so 
~ (O — q;)"'^/® in the critical region. The appropriate boundary condition 
is then V{0,r]i,Y) = 0. 

The propagator V{ri2,rii,Y) can be expressed by normalized eigenfunctions 
ipniv) eigenvalues En of the equation 

Eni^n = + {3v/2f/'^n • (3.10) 

Then 

V{v2,Vi,y) = J2Mm)MV2)e-^"'' . (3.11) 

n 

The probability distribution of the mid-point height is given by ipoiv) 
other quantities can be expressed by the appropriate combinations of eigen- 
functions. 

IV. The short-distance correlation function 

Obviously the above one-dimensional approximation cannot describe the 
full two-dimensional structure of the interface. However, there are two dif- 
ferent length-scales in this problem. The one- dimensional character of the 
filling transition is seen on scales ay ~ T,i'^ while the two-dimensional struc- 
ture becomes important when ay ~ i. In the critical region these two scales 
are well separated. 

Therefore, in order to analyze the short- distance behavior one can introduce 
the conditional correlation function. This is done in the following way. We 
assume that for certain yo (for convenience we set yo = 0) the interface pro- 
file i{x,y) is constrained: i{x,yo = 0) = i{x), where i{x) is described by 
Eq.( p.4|) but with given i^. The full Hamiltonian is then expanded in the 
Taylor series in the variable (f){x, y) = £{x, y) —i{x, 0) up to 0^ terms. In this 
way one obtains (up to the constant term) 



H[<P] = -2 dy 



aS - V2S(u;(4) - cu(£oo)) 0(0,1/) 

+ -jdxjdy [E(V0)2 + a;"(£»)02]. 

The first term on the rhs of Eq.(4.1) is very small in the critical region (i.e. 
for Iq ^ CO and 9 ~ a); it is given by ~ 211(6 — a). Thus for short 
distances one keeps only the second term. The resulting structure of the 
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Hamiltonian implies the following differential equation for the conditional 
correlation function 



G,„(r,r') = (0(r)0(r'))|,(.,o)=^c 



x) 



[-EAr + u"{e)] G,„(r, r') = 6{r - r') . (4.2) 

Similarly as in Eq.(3.11) the conditional correlation function can be expressed 
by the normalized eigenfunctions ip^ and eigenvalues of the operator 
[-EA + u;"( 



G,„(r,r') = ^ 



E, 



q 



In this approach one has to analyze carefully the contribution from the eigen- 
values tending to 0. One expects that the eigenfunctions with the lowest 
eigenvalues will have their structure similar to ipo = which itself corre- 

sponds to the translational mode (although it does not satisfy the appropriate 
boundary condition in the present case) . Thus we introduce the new variables 
■0q = (fiqipo and the equation for (pq has the form 



[Eq + EA] 



^2E(a;(£)-^(£oo)) 



d^ipq . (4.3) 



The expression on the rhs of the above equation is close to for Q;|a;| < £o and 
for a\x\ approaching £o it quickly becomes equal to —2^yTIU}"{£oo)^x'Pcl■ We 
are interested only in the long- wave fluctuations such that Eq ~ T,{a/£o)^. 
If all terms in the above equation are to be of the same order of magnitude 
for a\x\ > then one should have dx\Tnpci ~ Eq^T./T. ~ E,TvCi^ / (-1, where 
= ('-^"(^oo)/^)"^''^ is the correlation length for the planar case. Note that 
which diverges at the critical wetting on the planar substrate remains fi- 
nite at the critical filling transition. 

The above considerations lead to the following equation for the conditional 
correlation function for a\x\ < £o, i.e for the central "free" part of the inter- 
face 

-EArG,„(r,rO = S{r2-r,) (4.4) 
9,G',„(r,r')||.|=Va = 0. (4.5) 

One also needs the boundary condition for y' ^ oo. Since there exists the 
long-range order on the scale considered now, i.e. for ay ~ £o one should not 
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expect (^^^(ri, ^^^^ Q, Instead we assume (j£o(ri, ^^^^ /(ri) < oo, i.e. 
G remains finite. Using the standard methods of conformal transformations 
(see Appendix) one obtains the following solution of Eqs.(4.4, 4.5) 

G,o(ri,r2) = --^ [In (e^^-^^^)-/^ + e^^-^^)'^/^ - 2cos(Xi - X,)n/2) 

+ In (e(^i-^^)-/2 + e(^^-^^)-/2 + 2 cos(Xi + X2)7r/2) 

- In (e(^i+^2)^ + 1 - 2e(^^+^^)^/2 ^^^^^i - X2)7r/2) 

- In (e(^i+^2)" + 1 + 2e(^i+^2)"/2 ^^^^Xi + X2)7r/2) + n{Yi + Y2)] 

(4.6) 

where Rj = ari/io, i = 1, 2. We note that for r2 00 one has Gig{ri, r2) 

V. The short-distance dispersion 

For short distances the two-point ^-distribution function has the form 

p(^.. r, r.) « p(^.)!^E(z|i^ill!ZM£^ (5.1) 

[27rcr(ri,r2,£o)J^/^ 

where ^Q is the height of the interface above the edge of the wedge 
^0 = (^1 + 4 + "(kil + |a;2|))/2 ^ £1 + ^ £2 + a|a;2| ■ 

We use the conditional correlation function to obtain the expression for 
the dispersion a 

o = m-iif) = G,„(ri,ri) -2G,„(ri,r2) + G,„(r2,r2) . (5.2) 

The standard problem which one encounters at this point is that Gig{ri, r2) 
diverges for r2 ri [jl5|, |16| . This divergence can be removed by regularizing 
the function Ge^i^n, ^2), e.g. by adding to the Hamiltonian given in Eg. ([4.1|) 
the term a^(A0)^/2, where a is a dimensionless parameter. This procedure 
yields the following equation for the regularized function G['^\ri, r2) 

[a'Al^ - SArJ G[:^ = 6{r, - r,) . 
For small a the solution of the above equation has the form 

Gj;;)(ri,r2) = G,„(ri,r2) - MJ^'^'lr^ - ri|/a)/27rS , 
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Figure 2: The dimensionless dispersion a as function of Xi and F for X2 = 
and 0.5, respectively. 

where Kq is the modified Bessel function. In this way the short-distance 
divergence is removed and one has 



where 7 is the Euler constant. Now the expression for the dispersion a(ri, r2, ^o); 
Eq.(5.2) can be written down exphcitly. We are interested in the situation in 
which the constraint affects only the mean height of the interface and thus 
we consider the case yi,y2 ^ ^o/a and |ri — ^ a. Then 



(T(ri, r2, l^) = -—{2\n / aTia) + 27 - In (cos(Xi7r/2))) 

- In (cos(X27r/2)) + In [ch((Fi - ¥2)11/2) - cos((Xi - X2)7r/2)] (^-3) 
+ In [ch((ri - r2)7r/2) + cos((Xi + X2)7r/2)]} . 



The behavior of a is shown on Fig.§. We see that in this limit (j(ri,r2,^o) 
depends - in addition to Xi and X2 - only on the distance y = 1^2 ~ ^i- For 
fixed values of Xi and X2 it is an increasing function of \Y\, see Fig. 2. Thus 
the relative fluctuations of the interface position at points distant along the 
edge of the wedge become large. 

It is interesting to observe that for — ?/2| ^ ^o/a one gets 



This result agrees with the prediction of the one-dimensional model valid on 
the scale where S£g ^ — ?/2|- It can be derived with the help of Eq.(3.9). 



G'J°^(ri,ri) = lim G'^o(ri,r2) + 

" r2-»ri 



7 + ln(SV^|r2-ri|/2a) 
27rS 



a(ri,r2,4) ~ bi - y2|a(S£o)"V2 • 
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Thus the results obtained via the conditional correlation function in Chap- 
ters IV and V are consistent with those stemming from the transfer-matrix 
analysis of the ID model in Chapter III. 

VI. Conclusions 

The reduced description of the interface fluctuating in the presence of 
the wedge-shaped substrate is derived in an explicit way. This reduced de- 
scription is based on the one-dimensional Hamiltonian |jTO| and the presented 
derivation of this Hamiltonian makes clear use of the physical assumptions 
behind it. Although the one-dimensional Hamiltonian allows one to find 
the relevant critical exponents it cannot describe the full two-dimensional 
structure of the interface. We have proposed the method of supplementing 
this one-dimensional picture by the local two-dimensional constrained fluc- 
tuations which can be analyzed within the mean-field theory and described 
by the conditional correlation function. These fluctuations are not divergent 
at the filling transition. The proposed method can be used to calculate the 
geometry-dependent observables. Moreover, it predicts the behavior of the 
dispersion of the conditional correlation function which agrees with the pre- 
dictions of the one-dimensional model. 
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Appendix 

In this Appendix we sketch the consecutive steps of the method of con- 
formal transformations which lead to the solution of Eqs. (4.4,4.5). 
After introducing the complex variables ri^2 = (a;i,2, ^1,2), -21,2 = xi^2 + iyi,2 
Eq.(4.4) can be rewritten as 

-4S9,,9,,G,„(^i,^2) = 5(^1 - Z2). (A.l) 

together with the boundary condition (Eq.(4.5)) 

i[dzcd^ - dzcdz] Ge^iz, ^2) = , (A.2) 

where Zc denotes the contour on which the boundary condition is given. 
Eqs. (A.l, A.2) are invariant with respect to the conformal transformations. 
The solution of Eq.(A.l) valid for the whole plane has the form 

Goc{zi, Z2) = --^[log(2;i - Z2) + log(zi - Z2)]. (A.3) 
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The solution valid for the semi-plane > with the Neumann condition 
on^z — is found with the help of the method of images and has the form 

Gs}i^>o{zi, Z2) = Goo(^i, Z2) + Goo{zi, -Z2). (A.4) 

After introducing the conformal transformation Z 1— > e~^^^^^ for dimension- 
less variables Z — za/£ ,11— ra/£ one obtains 

Z2) = GKz>o(e^"^^/^ (A.5) 

with the Neumann condition at X = ±1. Finally, taking into account the 
constraint imposed on the interface at y = and using the freedom to add to 
the rhs of Eq.(A.4) the solutions of the Laplace equation leads to the solution 
of Eq.(4.4) given in Eq.(4.6). 
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